Nuclear spin relaxation in ordered bimetallic chain compounds 
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A theoretical interpretation is given to recent proton spin relaxation-time (Ti) measurements on 
NiCu(C7H6N206)(H20)3-2H20, which is an ideal one-dimensional ferrimagnetic Heisenberg model 
system of alternating spins 1 and |. The relaxation rate Tj -1 is formulated in temrs of the spin- 
wave theory and is evaluated by the use of a quantum Monte Carlo method. Calculations of the 
temperature and applied-field (H) dependences of T-j -1 are in total agreement with the experimental 
findings. Ti behaves as Tf 1 oc H 1//2 , which turns out an indirect observation of the quadratic 
dispersion relations dominating the low-energy physics of quantum ferrimagnets. 

PACS numbers: 75.10.Jm, 75.40. Mg, 75.30. Ds, 76.60.Es 



Recently considerable attention has been directed to 
ferrimagnetic mixed-spin chains. The simplest but prac- 
tical example of a quantum fcrrimagnct is two kinds 
of spins S and s alternating on a ring with antiferro- 
magnetic exchange coupling between nearest neighbors. 
Because of the noncompensating sublattice magnetiza- 
tions, this system exhibits ferrimagnetically degenerate 
ground states |Q. Thus, in contrast with ferromagnets 
and antiferromagnets, ferrimagnets show the ground- 
state excitations of dual aspect j|. The elementary ex- 
citations of ferromagnetic features, reducing the ground- 
state magnetization, form a gapless dispersion relation, 
whereas those of antiferromagnetic features, enhancing 
the ground-state magnetization, are gapped from the 
ground state. The two distinct low-lying excitations 
result in the unique thermal behavior which contains 
both ferromagnetic and antiferromagnetic aspects 
Quantum ferrimagnets in a magnetic field provide further 
interesting topics such as the double-peak structure of the 
specific heat |?J and quantized plateaux in the ground- 
state magnetization curves i 10 . 



The first ferrimagnetic one-dimensional compound, 
MnCu(dto) 2 (H 2 0)3-4.5H20 (dto = dithiooxalato = 
S2C2O2), was synthesized by Gleizes and Verda- 
guer 10]. Their pioneering efforts, combined with 
semiclassical but enlightening calculations, success- 
fully characterized the one-dimensional ferrimagnetic 
behavior fl2]| . Another example of an ordered 
bimetallic chain, MnCu(pba)(H 2 0) 3 -2H 2 (pba = 
1, 3-propylenebis(oxamato) = C7H6N2O6), was provided 
by Pei et al. |i~3| and turned out to exhibit more pro- 
nounced one dimensionality. Kahn et al. fli~4| further syn- 
thesized the family of related compounds systematically, 
focusing on the problem of the crystal engineering of a 
molecule-based ferromagnet— the assembly of the highly 
magnetic molecular entities within the crystal lattice in a 
ferromagnetic fashion. There also appeared an idea |l5| ] 
that the alternating magnetic centers do not need to be 
metal ions but may be organic radicals. 

In comparison with the accumulated chemical knowl- 
edge and well-revealed static properties introduced 
above, the dynamic properties of quantum ferrimagnets 



have much less been studied so far. To the best of our 
knowledge, it was not until quite recently that the dy- 
namic structure factors were calculated Jig ] , whereas any 
direct observation of the dispersion curves is not yet so 
successful, for instance, as that for the Haldane anti- 
ferromagnets |p^| . While Caneschi et al. jlfj performed 
electron-paramagnetic-resonance measurements on novel 
ferrimagnetic compounds, which consist of alternating 
metal ions and stable organic radicals, they focused their 
analyses on the relation between the spectra and short- 
range correlation effects. In such circumstances, the 
proton spin relaxation time T\ has quite recently been 
measured [[l9| for the spin-(l, i) ferrimagnetic Heisen- 
berg chain compound, NiCu(pba)(H20)3-2H 2 0, where 
the dependence of the relaxation rate on the ap- 
plied magnetic field H looks like T-j - oc H~ x /' 2 , though 
the authors' semiquantitative approximate argument has 
reached a logarithmic behavior, T^ 1 ~ \nH, rather than 
the observations. Motivated by this stimulative experi- 
ment, we here discuss the nuclear spin relaxation pecu- 
liar to ordered bimetallic chain compounds in connection 
with their energy structure. Formulating the relaxation 
process based on the spin-wave theory, we perform nu- 
merical evaluation of T^ 1 . Investigating the dependences 
of T x on temperature and the applied field, we demon- 
strate that the T\ measurement is nothing but an indirect 
observation of the characteristic elementary excitations 
of quantum ferrimagnets. 

The ferrimagnetic material NiCu(pba)(H20)3-2H20 
|]i~3f consists of ordered bimetallic chains with alternating 
octahedral Ni 2+ and square-pyramidal Cu 2+ ions bridged 
by oxamato groups. The chain runs along the 6-axis of 
the orthorhombic lattice whose space group is Puma and 
is described by the Hamiltonian 
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JV 



(1) 

where we have set the g factors of spins 5=1 and s = 5 
both equal to g because the difference between them 
amounts to at most several per cent of themselves in prac- 
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ticc pCfl . The high-temperature susceptibility measure- 
ments suggest that J ~ 121 [K]. The one dimensionality 
is well exhibited above about 7[K]. Due to the energy- 
conservation requirement for the electronic-nuclear spin 
system, the direct (single-magnon) process is of little sig- 
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nificance but the Raman (two-magnon) process plays a 
leading role in the nuclear spin-lattice relaxation 
Neglecting the higher-order relaxation process, the re- 
laxation rate is generally represented as 



Y, t iA'i S i + a i S i)l n )l 6 ( E m -E n - hu, N ) , 
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where A z and a z are the dipolar coupling constants be- where = (a\.ak) and n k = {fi\(3k) are the ther- 
mal averages of the numbers of the spin waves at a 
given temperature, and A z q = e 19< - J_1 / 4 ' and a z = 



tween the proton and electron spins in the jth. unit cell, 
cjn = 7N-ff is the Larmor frequency of 1 H with 7n be- 
ing the gyromagnetic ratio, and the summation J^n * s 
taken over all the eigenstates of \n) with energy E n . 
Let us introduce the bosonic spin-deviation operators 



via SJ = 
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Sj = — s + b'bj. The Bogoliubov transformation 

a k = cosh6> fc afe + sinh0 fc b' k , 
f3 k = sinhflfc a\ + cosh6» fc b k , 

with 
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diagonalizes the Hamiltonian (|l|) , discarding the constant 
terms, as 0J2T 



(6) 



where a\ and /?| are the creation operators of the ferro- 
magnetic and antiferromagnetic spin waves of momentum 
k whose dispersion relations are given by 



wjj = J\j (£- s) 2 + 4S*ssin 2 (fc/2) 
± (S-s)J T9^bH, 



(J) 



with twice the lattice constant being taken as unity. In 
terms of the spin waves, the relaxation rate (||) is ex- 
pressed as 
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it/(i+i/4) a z are ^ ne F our i er components of the hyper- 



fine coupling constants. Taking into account the signifi- 
cant difference between the electronic and nuclear energy 
scales (hujN < 10~ 5 J), Eq. d|) ends in 
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x [(A^cosh 2 ^. - e^sinh 2 ^) 2 ^ (n k + 1) 
+ {A z smh 2 k - a z cosh 2 e k ) 2 n+(n+ + 1)] , 



(9) 



where we have assumed little fc-dependence of A z and a z , 

A z =0 and 



and thus replaced A z = _ 



a' 



a„ 



respectively. 
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Now we have an idea of evaluating the thermal aver- 
ages rvt by a quantum Monte Carlo method. Though 
the divergence of the ground-state sublattice magneti- 
zation, which plagues the spin-wave treatment of low- 
dimensional antiferromagnets, does not persist in our 
model, n k still diverge as temperature increases as far as 
we naively estimate them with the noncompact bosonic 
Hamiltonian (|^). On the other hand, we are fully allowed 
to rely upon the dispersion relations (Q) themselves even 
at finite temperatures, because we can efficiently modify 
the spin-wave thermodynamics 0] so as to describe the 
higher-temperature behavior introducing an additional 
constraint on the magnetization but preserving the lin- 
earized spin-wave dispersions Wt. Thus, depending on 
the Bogoliubov transformation (|3|) but avoiding the di- 
rect estimation of the boson numbers, we perform the 
Monte Carlo sampling for the relevant spin operators 
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with the original compact Hamiltonian ([!]). While we 
base the relaxation process on the spin- wave excitations, 
we take grand-canonical averages within the original sys- 
tem. Since the applied field H is so small as to sat- 
isfy g/iB-ff ^ 10~ 2 J, we neglect the Zeeman term of the 
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in mono-spin chains, now shifts to the zone center due to 
the double periodicity in the present system. With the 
increase of temperature, the structure factor generally 
decreases in the vicinity of its peaks, while it increases in 
its lower slopes. Thus, due to the predominance of the 
k = component in the summation (g), T^ 1 ends up as a 
decreasing function of temperature. The increase of the 
applied field, having an effect of reducing the predomi- 
nance of the k = component, relaxes the pronounced 
decreasing behavior of T^ 1 as a function of temperature 
and would possibly allow T^ 1 to show a reversed but 
weak temperature dependence at high temperatures. 

Another interesting observation is the field dependence 
of the relaxation rate shown in Fig. |j. T^ 1 looks al- 
most linear with respect to _ff -1 / 2 , which also suggests 
that only the magnons of k ~ effectively contribute to 
the relaxation process. We stress that the present field 
dependence comes from the energy-conservation require- 
ment S(E m — E n — fno^) in Eq. (Q), that is, from the 
quadratic dispersion relations of the elementary excita- 
tions. In this context, the spin-diffusion effect |2j| may 
be mentioned, which also results in a similar field depen- 
dence of T-T 1 , 
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H 
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where the constant P arises from longitudinal spin fluctu- 
ations, whereas the field-dependent term originates from 
transverse spin fluctuations which are diffusive at high 
temperatures and long times. For (CHa^NMnChj |2J] 
and L1V2O5 p5| , which are both well known to be spin- 
diffusive Heiscnberg chain compounds, P/Q was esti- 
mated to be 0.57[T~ 1/2 ] at 300[K] with J = 13[K] and 
1.53[T~ 1/2 ] at 200[K] with J = 308[K], respectively. Here 
we observe a much smaller value, P/Q = 0.05[T~ 1//2 ] at 
280 [K], which is much more consistent with the present 
analysis. Therefore, the present observations, which look 
like T-f 1 cx H^ 1 / 2 , should be attributed to the quadratic 
dispersion relations peculiar to Heisenberg ferrimagnets. 
We learn that the longitudinal spin fluctuations dom- 
inate the nuclear spin relaxation in ordered bimetallic 
chain compounds. 



The present compound, NiCu(pba)(H20)3-2H20, is 
reasonably described by the isotropic Heisenberg Hamil- 
tonian. There exist a series of family compounds Jll] 14] 
and thus the prototypical ferrimagnetic behavior 26] 
could potentially be elucidated as a function of the con- 
stituent spins. The consideration of anisotropy effects 
may also be interesting from the practical point of view. 
For the anisotropic Hamiltonian 
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with (S ■ s) a = S x s x + S y s v + aS z s z , the relaxation rate (j|) is replaced by 
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where 



f s , s (a,D s ,D s ) = y/[a{S + s) - (SD S + sD s )/J] 2 - ASs . 
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Even if any anisotropy of the easy-axis type is introduced, 
the small-momentum quadratic dispersion remains qual- 
itatively unchanged and thus the characteristic field de- 
pendence of the relaxation rate, Tj~ cx H^ 1 / 2 , should 
still be observed, where the slope dT^ 1 / dH^ 1 / 2 gives 
semiquantitative information on the anisotropy. How- 
ever, the expression ( |l3| ) is no more valid under the easy- 
plane-type anisotropy. The model indeed turns critical 
for a < 1 , showing linear dispersion relations j p7| . There- 
fore ferrimagnetic compounds with anisotropic interac- 
tions of the easy-plane type should exhibit no field depen- 
dence of T±. Thus the nuclear magnetic resonance is all 
the more efficient to investigate the low-energy structure 
in the present system. We really hope that the present 
successful interpretation of the pioneering experiments 
will stimulate and accelerate further study on this fasci- 
nating system. The developed calculations under various 



constituent spins and geometric parameters will be pre- 
sented elsewhere p8[ . 
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